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Abstract. 

Quantum fluctuations in the mazer are considered, arising either from the atomic motion or from the 
quantized intracavity field. Analytical results, for both the mcza and the hyperbolic secant mode profile, 
predict for example an attenuation of tunneling resonances due to such fluctuations. The case of a Gaussian 
mode profile is studied numerically using a wave packet propagation approach. The method automatically takes 
into account fluctuations in the atomic motion and the dynamics is especially considered at or adjacent to a 
tunnel resonance. We find that the system evolution is greatly sensitive to the atom-field detuning, bringing 
about a discussion about the concept of adiabaticity in this model. Further, a novel collapse-revival phenomena 
is demonstrated, originating from the quantum fluctuations in the atomic motion rather from field fluctuations 
as is normally the case. 



PACS numbers: 42.50.Pq, 03.65.Xp, 32.80.Qk, 03.75.-b 

1. Introduction 

In recent years, the field of experimental cavity quantum electrodynamics (QED) has expanded beyond the 
setup of single thermal two-level atoms traversing the cavity one by one [U [2 |3j . Cavity QED now includes 
the realization of solid state q-dots coupled to cavity fields [4], trapped ion-cavity systems [5], Bose-Einsten 
condensates coherently interacting with single cavity modes \§\ and ultracold atomic many-body systems inside 
cavities; e.g., bistability behavior [7], cavity cooling [5j and self-organization of atoms in regular patterns [S]. 
Reaching these new regimes of cavity QED has been possible due to the progress in the cooling, trapping and 
coherent control of atoms [TO]- Actually, atomic velocities of the order of 1 cm/s or even slower are achievable 
[11] . This opens up the possibility to experimentally demonastrate the mazer (Microwave Amplification via 
z-motion-induced Emission of Radiation) , where ultracold, rather than thermal, two- level atoms are let to pass 
the cavity [121 [13]. 

In the ultracold regime, the kinetic energy of the atom is of the order of the atom-field interaction energy 
and must be treated quantum mechanically. Early on, it was predicted that atoms can be back scattered by 
the high-Q cavity vacuum induced potential [12j or even trapped by the intracavity field 14], which indeed was 
later experimentally obtained for single photon fields [15 . Later, in a series of papers by the late Walther and 
coworkers, it was shown that in this new regime the emission probability for an excited atom traversing the cavity 
is greatly affected by the quantized atomic motion [13] . For certain resonance conditions, where the atomic de 
Broglie wavelength match resonantly with the cavity length, the atom tunneling and emission probabilities are 
greatly increased. Such resonances lead to new phenomena of the field statistics in the micromaser pumped by 
excited atoms. Following [13] . a set of papers where published on the mazer. Most of which consider analytically 
solvable mazer models, namely the atom-field detuning vanishes and the cavity mode profile has the shape of 
a meza function or a hyperbolic secant [TBI HZl [18 , 19, 20 . These all assume the zero detuning situation where 
the mazer problem can be decoupled into two separate 1-D scattering problems, while for a non-zero detuning 
the dynamics is indeed more complex and the full coupled system must be considered. This was the subject of 
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[21] , where the problem of a non-zero detuning in the case of a meza function coupling was solved analytically 
(see also [22 ). Further, the situation of a sinusoidal mode profile with one or a few wavelengths and a Gaussian 
shape have been analyzed by means of approximate or numerical methods in [23] and [18] respectively. The 
limit in which the mode profile is "delta" -shaped has as well been discussed [231 [33]. The mazer problem for a 
sinusoidal coupling with several wavelengths was investigated in [26] and the model was found to exhibit very 
different characteristics due to the quasi-periodicity, see also [27] . Extensions of the regular mazer setup have 
been considered in numerous works, e.g. a driven system [28] . multi-level atoms and multi-photon processes 
[29] , multi-cavity setups [El [22J, [30] , Kerr and gravity effects on the mazer action j3TJ [32] . Also semi-classical 
approaches have been studied, where for example the atomic motion is replaced by an effective time-dependence 



None of the above references take into account fluctuations in the motion of the atom. In the standard 
atom cavity QED experiment in which the atom traverses the cavity, its position is fairly well localized in space, 
typically for a microwave setup within the order of 1 mm (where the mode waist L is about 5-6 mm) pQ. Thus, 
in real experiments is the state of the atom given by a well localized wave packet rather than a traveling wave. 
Naturally, an uncertainty in position results in fluctuations in momentum which, on the other hand, affects the 
particle scattering. In this paper we address these issues both by generalizing known analytical models but also 
via numerical analysis of a Gaussian mode profile. The Gaussian coupling approximates the mode profile of a 
Fabry-Perot cavity which is usually used in micromaser experiments [TJ E] . A crucial feature of the Gaussian, 
as well for the hyperbolic secant, coupling is that they are smooth with a well defined analytical limit, contrary 
to the meza function. This give rise to qualitative different results compared with the meza case. Here we 
focus on these consequences and discuss the concept of adiabaticity in the mazer. One essential parameter for 
adiabaticity is the atom-field detuning and how it affects the dynamics is investigated. It is found that it plays 
a central role for the evolution. So far, non-zero detuning has only been considered for meza-mode couplings, 
which do not possess an adiabatic limit due to their discontinuities. Quantum fluctuations lead in general to a 
smearing out of physical observables, but we demonstrate that it can also render new collapse-revival structures. 

The general mazer concept and the model system are introduced in the section[2l The various reflection and 
transmission coefficients, as well as appropriate bases, are given. Our analytical results are presented in section 
[3j where quantum fluctuations in both the meza and hyperbolic secant cases are considered. The following 
section [4] is devoted to numerical research of the Gaussian mode profile for zero and non-zero detunings. We 
apply a split-operator approach that invitably takes into account fluctuations in the motion of the atom. The 
effect of the detuning parameter is studied as well as the system dynamics during the scattering process. In 
the final section [5] we give a summery and discuss experimental prospects. 

2. The model system 

We consider a two-level atom with internal states \g) and |e) and center-of-mass (c.o.m.) momentum and 
position p and z, which dipole interacts with a single quantized mode of a cavity. The effective Hamiltonian in 
the rotating wave approximation reads 



where m is the mass of the atom, u> the field frequency, f2 is the atomic transition frequency and X(z) is the 
effective position dependent atom-field coupling. Here a and are the annihilation and creation operators for 
the field and the u-operators are the regular Pauli matrices, a z = |e)(e| — cr~ = \g){e\ and a + — |e)(<?|. 

The number of excitations is conserved and it may be used to define an interaction picture rotating with uj. This 
introduces the atom-field detuning A = £1 — uj and the transformed Hamiltonian becomes (we let K = m = 1) 



[31 [34]. 




(1) 




(2) 
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Due to this symmetry can the Hamiltonian be written on 2 x 2 block diagonal form within its internal states. 
For a given number of excitations, corresponding to one of the 2 x 2-blocks, it is often convenient to turn to a 
dressed representation by transforming H with [35] 

cos(0) sin(0) 1 , s 

[ - sin(0) cos(0) J ' [ ' 

where tan(2#) = 2A(z)y / n/A. This defines the dressed internal states 

| 7 +> = cos(0)|n,e> + sin(0)|n + l,s), 

(4) 

|7~) = - an(6)\n, e) + cos(9)\n + 1, g). 

Here, \n, e) and \n + 1, g) are the states with n photons and the atom excited or n + 1 photons and the atom 
in its ground state respectively. These states are referred to as bare internal states. The bare state |0,g) is 
decoupled from any other states, which, however is not true beyond the rotating wave approximation. We 
further have that U diagonalizes the last two terms in the Hamiltonian ([2]) whose eigenvalues read 




n = 0,l,2,.... (5) 

We call these adiabatic potentials whenever A ^ 0. Since the operator U is z-dependent it will not commute with 
the first term of H, causing non-diagonal terms which couples the two dressed internal states [35]. However, 
in the limiting cases A — > ±00, one has that these non-adiabatic couplings vanish. Thus, for A = we may 
decouple the two equations and we note that in this special case cos(#) = sin(#) = 1/V2. We distinguish this 
limit with the ones of large detuning (adiabatic limit) and call it diabolic limit [36] . 

Due to the coupled two-level structure, the dynamics may become rather complicated. An atom entering 
the interaction region will in general feel both an attractive and a repulsive potential. Despite the coupled 
dynamics, transmission and reflection coefficients for the two internal states can be well defined. In particular, 
for the dressed basis we call these p^ and ■ It follows that the coefficients in the bare basis are related to 
the dressed ones via [IJJ 

R en = ~(pt+Pn), R 9 n = \{pt~ Pn ) > 

I l" 

T en = 2 { T n + T n ) ) ^gn — ^ { T n — T n ) ■ 

The i?-coefficients squared give the probability for reflection of the atom in state e or g respectively and 
correspondingly T gives the transmission. Thus, we have the transmission probability, given an excitation n, 

Pfrans = \T e n\ + \T gn \ (7) 

and similarly for the reflection. It should be pointed out that the mazer scattering process is different from 
the regular 1-D ones often presented in textbooks [37]. This comes about because of the coupled dynamics 
between the internal states, and in particular, in some limiting cases can the particle be said to be scattered 
simultaneously to an attractive and repulsive potential. This, in fact, leads to new phenomena [16] . not 
encountered in the standard situations. Another property of interest is the von Neumann entropy which gives 
a measure of entanglement shared between the field and the atom. Defining the reduced density operator for 
the field as pf = Tr Qt [/o], where p is the full system density operator and the trace is over atomic degrees of 
freedom, one has the von Neumann entropy 

S n = -Tr[p f ]0g(j> f )]. (8) 
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The field and atomic entropies are identical for initial pure states [35 
that 

S n = -(|i? CTl | 2 + |T en | 2 )log[(|i? en | 2 + |T e „| 2 )] 



(|i? gn | 2 + |T s „| 2 )log[(|i? 9 ' 2 



In terms of the coefficients © it follows 



(9) 



. . </" i " U " '/" ! r |7g„| 2 )] . 

Let us choose the initial state to be the uncorrelated state with the field in vacuum, the atom in its internal 
excited state and with some spatial wave function ip(z), 

(*|tf(t = 0)) = Vi(*)|0,e>. (10) 
Using the position-momentum relation 

i>(z) = J dktP(k)e lkz (11) 

we can write down the wave function for large times, when the atom is far from the interaction regime, as |13j 

(*|*(t)) - [ dk^e-^/^ilReoi^e-^ +T g0 (k)e ikz ]\0e) 



(12) 



[R g0 (k)e- ik *+T g0 (k)e ik *] \l,g)}. 



Taken into account for fluctuations in momentum we define the total reflection and transmission coefficients 



Rio 



dk 



dk 



^(k)R lQ {k) 



ip(k)T i0 (k) 



e, g, 



e, g. 



(13) 



With these, the transmission probability or the von Neumann entropy are readily modified. Noteworthy is that 
the coefficients (fTB"]) is as well easily generalized for general initial cavity fields, \<p) = J2 n Cn \ n ) by weighting 
the corresponding coefficients by \c n \ 2 and summing over n. 

3. Analytical consideration 

Some of the outcomes from quantum fluctuations, both in the atomic c.o.m. motional state and in the cavity 
photon field, are demonstrated in this section using analytically solvable models. The sequential section 
continues the analysis by a numerical study of a more physically realistic situation. In both models considered 
here is the detuning A is assumed zero. In other words, the internal structure can be separated into two 
decoupled equations. The atomic motional wave function is taken to be a Gaussian, centered at kg and with 
width Afc, 



1 



-A 2 



1/4 (k-kp) 2 
„ 4Af 

e fe 



(14) 



If not mentioned, the cavity field is initially in the vacuum. 



3.1. Meza function 

The simplest solvable model, and also most extensively studied, assumes a constant atom-field coupling within 
and zero outside the cavity region p~2l [13j [TBI ttZ] 



X(z) = 



A 




for < z < I 
elsewhere. 



(15) 
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Figure 1. The transmission probability Ptrans and the field entropy S (for the meza function) as function of 
atomic c.o.m. momentum fco- Solid lines correspond to the results J7J and while dashed and dotted lines 
take into account for momentum fluctuations. Here Ao = 1 and L = 50. 



The dressed transmission and reflection coefficients read [13] 

Pn = iA t sin(>± L)r±, 

(16) 

T t = expHfci) [cos(fc±L) - jE± sin(fc±L)] , 

where 




K n = Ky/n~+ 1 . 

The discontinuity of the atom-field coupling, once the problem is presented in the dressed basis, will render 
delta-like non-adiabatic coupling terms [35 . This singularity naturally affects the system properties and in 
particular are such effects believed to be absent in realistic experimental models to be discussed in the next 
Section. 

The results for the transmission probability and the von Neumann entropy are presented in Figures [1] and 
[2j The first figure shows how the quantities depend on the atomic c.o.m. momentum kg. The solid lines 




give the zero fluctuation case (0 and Q, the dotted lines are calculated with a momentum fluctuation of 2 
% of the c.o.m. momentum, Afc = 0.02fco, and for the dashed line Afc = O.lfco- It is clear how the quantum 
fluctuations smear out any variations of the two quantities, including the tunneling resonances in the mazer 
regime < ko < \J~2. The atom-field coupling is here Ao = 1 such that for momentum ko > ko = V2Ao is 
the atomic kinetic energy larger than the potential barrier and the atomic transmission is abruptly increased 
beyond kq. In this Rabi regime one notes how the atom- field Rabi oscillations begin to form. The period of 
Rabi oscillations decreases faster than the width Afc increases, resulting in that these oscillations are captured 
in the dotted line (Afc = 0.02fco) for large k. The atom-field entanglement is maximized when the transition 
probability reaches 1/2 in the mazer regime. In the Rabi regime, the entanglement follows the swapping of the 
excitation between the field and the atom. In Figure [2] we show the same quantities when instead the coupling 
length L is varied. In this case, the results are insensitive to small momentum fluctuations and we therefore 
only show the results for Afc = 0.5fco- The resonances envisaged in the plots occur for L — mn/n n ^. Note 
that the fluctuations tend to increase the quantum correlations between the atom and the field. 

How an uncertainty in the field statistics affects the transmission probability is depicted in Figure [3J The 
upper plot (a) shows the transmission probability for an initial coherent field state 

|c„| 2 = ( 18 ) 
n! 

with mean photon number no = 0, no — 0-2 and uq — I. Figure 3 (b) gives the same but for an initial thermal 
field state 



| Cn | 2 = -^r— j-r. (19) 

1 1 (n + l) n+1 V ; 

It is found that the tunneling resonances are very sensitive to fluctuations in the cavity field. However, in 
the microwave regime of cavity QED, a thermal field with no = 1 (no = 0.2) correspond to a temperature of 
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Figure 3. The transmission probability Ptrans (for the meza function) as function of mode length L for an 
initial coherent cavity field (a) and an initial thermal cavity field (b). The average number of photons is no = 
(solid line), no = 0.2 or no = 1 (dashed line) and the other parameters are as in Figure|2] 



0.5 (0.2) degrees Kelvin, which is way above current effective experimental temperatures pQ. The secondary 
peaks showing up for non-zero no are the tunneling resonances arising from other photon quantum numbers n 
differing from n = 0. 

3.2. Hyperbolic secant function 

To simulate a more realistic situation we consider a smooth atom-field coupling 

A(z) = A sech 2 (|) , (20) 

where Ao is the effective atom-field coupling and L determines the mode waist length. The scattering problem 
analytically solvable in the special case of zero detuning A = 0, with dressed reflection and transmission 
coefficients [13] 



. r(i P L)r(i - ikL) ± 

r(i/2 + i£)r(i/2-i£±) 

(21) 



± = r[i/2 - ijikL + £)]r[i/2 - i{kL - g*)] 

n[ ' T(-ikL)T(l-ikL) 



where 



g - J±2X L 2 Vn+l - 1/4. (22) 



In terms of re n (as defined above), the tunneling resonances are given by n n L = ^Jm(m + 1) for m = 1, 2, 3, ... 
l.'ij. The transmission probability Ptt- ans (ko) and the von Neumann entropy S(ko) are displayed in Figure |4] as 
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Figure 4. The transmission probability Ptrans and the field entropy S (for the hyperbolic secant function) as 
function of the atomic c.o.m. momentum k,Q. Contrary to the meza function, the transition probability equals 
1/2 and the the atom-field becomes maximally entangled for most of the momenta ko in the mazer regime 
feo < Ao\/2. Here, as in Figure [T] Aq = 1 and L = 5. 



function of feo. The solid line gives the quantities for zero momentum fluctuations, while the dotted and dashed 
curves take into account for an uncertainty in the momentum distribution. It is seen that the transmission 
probability is constant in the mazer regime for this smooth coupling. In the Rabi regime is the transmission 
probability unity throughout composite to the meza function, which is an effect of the smoothness of the 
coupling. From the entanglement presented in plot (b) is it clear that in this regime the system Rabi oscillates. 
Comparing Figures [T] and [4J the effect of the non-analyticity of the meza function coupling is visible from the 
irregular behaviour of Figure [1] Again, the momentum fluctuations tend to blur the variations of the curves. 

The dependence of the cavity mode waist L is depicted in Figure [5] As L increases the transmission 
probability Ptrans(L) and von Neumann entropy S(L) approaches 1/2 and 1 respectively, which therefore 
deviate considerably from the results for the meza function presented in Figure O It is known that for large L, 
the adiabaticity is in general increased [35 . Here, however, do we consider the case of A = and we thus are 
in the strict diabatic limit rather than the adiabatic one. 

4. Numerical consideration 

In this section we consider a more realistic situation. The atoms perpendicularly traverse the field of a Fabry- 
Perot cavity, in which the mode shape is given by a Gaussian 1 

A(z) = A e"ir. (23) 

The case of a Gaussian atom-field coupling in the mazer was considered numerically in [IS], and the tunneling 
resonances where compared with the ones of the hyperbolic secant of the previous section. We, however, use 
a completely different numerical approach, namely a wave packet propagation procedure more commonly used 



Quantum fluctuations in the mazer 



9 



(a) i 
0.8 
fcl 0.6 




Figure 5. The transmission probability Ptrans and the field entropy S (for the hyperbolic secant function) as 
function of mode waist L. Contrary to the meza function, the tunneling resonances vanish for large L where 
the transition probability approaches 1/2 and the the atom-field system becomes maximally entangled, 8 = 1. 
Here, as in Figure [2] are Ao = 1 and fcrj = Arj/10. 

in molecular and chemical physics [39j . An initial wave packet, which in our case is taken to be a minimal 
uncertainty Gaussian 

= (d^) e _k ^i e '°>' ( 24 ) 

is propagated in time under the corresponding Hamiltonian (J2|) using the split operator method [JD]. Thus, the 
wave packet is obtained at any instant of time and physical quantities are easily extracted from it. In particular, 
the transmission probability is calculated as 

/>oo 

Ptrans = lml / dz \$( Z) t) | 2 , (25) 

where the time limit is taken such that the result has asymptotically converged. An advantage of this approach 
is that it automatically takes into account for motional quantum fluctuations via the width A z of the wave 
packet. Contrary to Reference [IS], we will be interested in the effect of having a non-zero atom- field detuning 
A as well as the full system dynamics during the scattering process. 

4-1. Tunneling 

As already pointed out, the meza mode profile is expected to give an incorrect descripsion of realistic 
experimental situations. Firstly, it does not have the actual shape of typical cavity mode profiles, and secondly, 
it is a non-analytic function. Indeed, it is known that discontinuities cause non-adiabatic behavior and cannot 
correctly reproduce realistic cavity QED situations [34J. More realistic results are instead expected to be 
obtained from the solvable hyperbolic secant model. However, this system is only solvable for a zero detuning 
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Figure 6. (Colour online) The transmission probability Ptrans and the von Neumann entropy S (for the 
Gaussian mode profile) as function of mode waist L and for three different detunings; A = (solid line), 
A = 0.02 (dashed blue line) and A = 0.1 (dotted red line). Wide mode profiles convey a more adiabatic 
dynamics, causing a decrease in Ptrans and S. Here, A z = 15, Ao = 1 and k = Ao/10. 

21j. Here we want to explore the scattering characteristics in the case of a Gaussian mode profile and compare 
with earlier results for the meza and hyperbolic secant. 

In Figure [6] we present the numerical results for the transmission probability Ptrans{L) and the von 
Neumann entropy S(L). The solid lines give the zero detuning case where the coupled equations can be 
decoupled into two scattering processes of attractive and repulsive potentials respectively. The tunneling 
resonances appears more often than for the hyperbolic secant coupling of Figure which was also pointed 
out in [IS]. Another difference is the more rapid convergence to the asymptotic values Ptrans(L) — 1/2 and 
S(L) = 1 for large L of the hyperbolic secant profile compared to the Gaussian. This is a consequence of the 
longer tails of the hyperbolic secant and the smoother behaviour in general. This feature has been studied in 
a semi-classical approach for cavity QED models in [33] , where it was found that the evolution is adiabatically 
more robust for a hyperbolic secant than for a Gaussian mode profile. However, it should be mentioned that 
the atomic c.o.m. motion was replaced by an effective time-dependence and the atomic kinetic energy exceeds 
the atom- field interaction energy in that work. The dotted and dashed lines in Figure [6] are the results obtained 
from having a non-zero detunings A. Note that in both cases A < Ao and still a great effect of the non- 
zero detuning is found. In other words, the tunneling dynamics is very sensitive to the parameter A, and we 
especially note that this detuning dependence is clearly stronger for a Gaussian mode profile compared to a 
meza function [5T]. We discussed the adiabatic and diabatic potentials in Section where |A| — > oo gives the 
adiabatic limit. Figure [6] confirms our earlier conclusions that an increase in adiabaticity depletes the tunneling 
effect; large A and L render an adiabatic evolution and in this regime the tunneling resonances are completely 
washed out. 
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4-2. Dynamics 

The atom-field dynamics in the mazer was briefly considered in |20j in terms of the characteristic tunneling 
time for a meza mode profile. Dynamics was also analyzed for an ultracold atom traversing a cavity with a 
cosine mode profile consisting of several periods [26] . It was found that due to the quasi periodicity an effective 
mass can be introduced in order to understand the evolution. The same system was investigated further in [25j 
studying the concept of a matter wave index of refraction. 




-600 -400 -200 200 400 600 



-600 -400 -200 200 400 600 



Figure 7. Time evolution of the probability distribution \ip e {z, t)\ 2 + \ijjg(z, t)\ 2 (a) and (b), and the individual 
amplitudes \if> e (z,t = 3000) | 2 (dotted) and \ip g (z,t = 3000) | 2 (solid) in (c) and (d) for a Gaussian mode profile. 
The left plots display the case where the parameters are chosen to give maximum tunneling (corresponding to 
the third resonance of Figure(6j while the right ones show the results with parameters in between two resonances 
(second minima of Figure[6]l. The insets zoom in on the wave packets in the classically forbidden region. Other 
parameters are as in Figure [6] 

We now turn to consider the Gaussian mode profile and the wave packet evolution while passing the 
interaction region, both when the tunneling condition is fulfilled and when it is not. We use the parameters of 
Figure [6] and chose L = 2.1245 or L = 1.8108 for the two situations respectively. These values of L correspond 
to the third tunneling resonance and the second minima of Figure El We restrict the analysis to the zero 
detuning case. The results from the wave packet propagation are depicted in Figure [71 First, in (a) and (b) we 
give the total wave packet amplitude \^(z, t)\ 2 = \ip e (z , t)\ 2 + \ijj g (z , t)\ 2 , while (c) and (d) show the constituent 
amplitudes \i/} e (z, t = 3000) | 2 (dotted) and \ip g (z, t = 3000) | 2 (solid) at time t = 3000. The insets give a close-up 
of the amplitudes around the interaction region, X(z) 0. The difference between the two examples of being 
in or ouside the tunneling regime is evident from the plots. The typical decrease of the wave packet inside 
the classically forbidden regime throughout the evolution is also visible from the figure [H]. The symmetric 
oscillations structure of the amplitude inside the forbidden region, seen in the inset of (c), is a characteristic of 
the tunneling resonance condition where the number of nodes depends on which resonance peak one considers; 
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the number of atomic de Broglie wavelengths that fit the cavity length. Note that some oscillations are also 
present in the forbidden regime also in (d), but one should bare in mind the difference in amplitude (factor of 
100) of the two insets. 

4-3. Motional induced collapse-revivals 

It is well known from quantum optics, and especially from cavity QED, that an uncertainty in the field state 
results in a collapse in the atom- field Rabi oscillations [HI 03] . The various photon number states that form 
the field state induce a series of different Rabi oscillations causing the collapse. At later times, the various 
terms comes back in phase, manifesting it selves in a revival. Also for the mazer has the collapse-revival effect 
been analyzed [44] . In our model do we consider an initial vacuum and no collapse- revival pattern can arise 
due to field fluctuations. We will, however, show that the fluctuations in the atomic motional state may form 
a collapse-revival structure. 




250 



t 

Figure 8. The atomic inversion W(t) and the von Neumann entropy S(t) (for the Gaussian mode profile) as 
function of time. Here A z = 10, Ao = 1 A = and fco = 6. 

Assuming k > K n , but remaining in the regime where k ~ K n , the whole wave packet will be transmitted 
through the vacuum induced potential. From Figured (b) we have that, while traversing the cavity field, the 
atom will Rabi oscillate with frequency A(z). However, since the coupling A(z) is not constant over the extent of 
the atomic wave packet, will the oscillation frequency vary within the packet. Entering the interaction region, 
the front of the wave packet sees a stronger amplitude of the atom-field coupling than the back of it. Letting 
L/2 and r = L/2kg denote characteristic length and time scales respectively, and further assuming that the 
wave packet is to the left of the origin 2 = and that A z < L. Then, typically if [A(0) — A(— A z )]r ~ it a 
collapse will occur in the dynamical quantities, such as the atomic inversion W(t) = (a z ) or the von Neumann 
entropy. After passing z = 0, the coupling decreases and as the whole wave packet is to the right of the centre 
the situation will be the opposite; the front of the packet oscillates with a smaller frequency than the back. 
Thus, because of the symmetry of the Gaussian, the collapse will be reversed resulting in a revival, similar to 
time-reversal echo processes [431 145j . The numerical results shown in Figure [5] for the inversion and the von 
Neumann entropy confirm the above arguments of a collapse-revival structer brought about due to momentum 
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fluctuations. The imperfect revival is because wave packet spreading during the passage of the interaction 
region. Note that the revivals occurs due to the "time-reversal" like dynamics and not because of a discreteness 
in quantum numbers, as is the origin of regular revivals [421 146j . 

5. Conclusion 

In this work we have investigated the mazer beyond earlier studies by allowing for quantum fluctuations in 
the atomic motion and also in the field statistics. Known analytical models, the meza or the hyperbolic 
secant, have been extended to include such uncertainties and their effects have been discussed. Typically the 
characteristics of various physical quantities, e.g. entanglement or transmission probability, are smeared out 
by these fluctuations. We have also shown that the fluctuations in general increase the amount of atom-field 
entanglement. 

The more physically interesting situation, in terms of experimental realization, of a Gaussian mode profile 
was considered via a numerical wave packet approach. This method naturally takes into account for motional 
fluctuations and it is not limited to the situation of zero atom-field detuning. A consequence of a non-zero 
detuning is the fact that a decoupling of the internal two-level structure is not straightforward, apart from the 
limiting cases of adiabatic or diabatic evolution. We found that the system dynamics is very sensitive to the 
value of the detuning, which opened up questions about adiabaticity contra tunneling behaviour. Using our 
numerical propagation procedure we also studied the dynamics of the wave packet during the scattering against 
the vacuum potential. Finally we demonstrated how motional fluctuation can induce novel collapse-revival 
characteristics, different from the regular cavity QED collapse-revivals originating from fluctuations in the field 
statistics. 

Experimentally, a difficulty with reaching the mazer regime is the extremely long characteristic time scales. 
The atoms are cooled down to the nK regime with typical c.o.m. velocity of mm/s |13j . Using microwave cavities 
of mm-size [11 implies operational times of a few seconds. This exceeds typical atomic and cavity life-times which 
are of the order of ms. Another obstacle for these long processes is gravity which however may be circumvented 
by using matter wave guides [47] . Atomic losses can be minimized by using a two-photon RAMAN coupling, 
such that two meta stable atomic states are dispersively coupled via an excited state using one external laser 
field and the cavity field [48]. This, however, introduces additional Stark shift terms in the Hamiltonian. To 
overcome the problem of cavity losses one could consider a pumped open cavity system. Another possibility is 
to try to speed up the process, for instance by considering a Bose-Einstein condensate coupled to the cavity 
[6], in which the effective atom-field coupling can be increased by a factor 1000 compared to the single atom 
case [49 . In such a system, the intrinsic non-linearity leads to novel phenomena like bistability. This is under 
current investigation. 
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